Effect of electron— phonon interaction on the shift and attenuation of optical phonons 
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Using the Boltzmann equation for electrons in metals, we show that the optical phonons soften 
and have a dispersion due to screening in agreement with the results reported recently [M. Reizer, 
Phys. Rev. B 61, 40 (2000)]. Additional phonon damping and frequency shift arise when the 
electron-phonon interaction is properly included. 
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Despite attracting considerable interest for half a cen- 
tury since the pioneering work by Frohlich, the prob- 
lem of electron-phonon interaction is still far from being 
solved. Migdal (l) developed a consistent many-body ap- 
proach based on the Frohlich Hamiltonian for interaction 
of electrons with acoustic (sound) phonons. As Migdal 
showed ("the Migdal theorem"), the vertex corrections 
for acoustic phonons are small by the adiabatic parameter 
y/m/M, where m and M are the electron and ion masses, 
respectively. The theory described correctly the elec- 
tronic lifetime, renormalization of the Fermi velocity vf 
and acoustic phonon attenuation but resulted in a strong 
renormalization of the sound velocity s = s(l — 2A) 1 / 2 , 
where A is the dimensionless coupling constant. For suf- 
ficiently strong electron-phonon coupling A — > 1/2, the 
phonon frequency approached to zero marking an insta- 
bility point of the system. Instead, one would intuitively 
expect the phonon renormalization to be weak along with 
the adiabatic parameter. 

This discrepancy was resolved by Brovman and Kagan 
M almost a decade later (see also ||). They demon- 
strated the shortcomings of the Frohlich model that gave 
an anomalously large phonon renormalization. Employ- 
ing the Born-Oppenheimer (adiabatic) approximation 
(see, e.g., Q), they found that there are two terms in 
the second order perturbation theory, which compensate 
each other making a result small by the adiabatic param- 
eter. Namely, when calculating the phonon self-energy 
function H(ui, k) with help of the diagram technique, one 
should eliminate an adiabatic contribution of the Frohlich 
model by subtracting U(cj, k) — n(0, k). 

The interaction of electrons with optical phonons 
was first considered by Engelsberg and Schrieffer || 
within Migdal's many-body approach for dispersionless 
phonons. They predicted a splitting of the optical 
phonon at finite wavenumbers k into two branches. Ipa- 
tova and Subashiev calculated later on the optical 
phonon attenuation in the collisionless limit and pointed 
out that the Brovman-Kagan renormalization should be 
carried out for optical phonons in order to obtain correct 
phonon renormalization. In the paper H, Alexandrov 
and Schrieffer corrected the calculational error of Ref. || 
and argued that no splitting was found in fact. Instead, 
they predicted an extremely strong dispersion of optical 



phonons, u>k = ujq + Xvpk 2 /3wo, due to the coupling to 
electrons. No such a dispersion has ever been observed 
experimentally. The usual dispersion of optical phonons 
in metals has the order of the sound velocity. In a recent 
paper, Reizer j| stressed the importance of screening ef- 
fects which should be taken into account. The works [Q, 
H are limited to the case of collisionless both electron 
and phonon systems. Moreover, only the phonon renor- 
malization was considered with no results available for 
the attenuation of optical phonons. 

A different from many-body technique semiclassical 
approach based on the Boltzmann equation and the equa- 
tions of the elasticity theory was developed in the papers 
by Akhiezer, Silin, Gurevich, Kontorovich, and many 
others (we refer the reader to the review |J). This ap- 
proach was compared with various experiments, such as 
attenuation of sound waves, effects of strong magnetic 
fields, crystal anisotropy, and sample surfaces on the 
sound attenuation, and so on. It can be applied to the 
problem of the electron-optical-phonon interaction 
as well. In the present paper we develop a theory for both 
the attenuation and frequency shift of optical phonons 
with account for effects of the Coulomb screening as well 
as collisions in the electron and phonon systems. 

It is instructive to recapitulate the results of the dy- 
namical theory of elasticity for the renormalization of 
the sound velocity 6s = s — s and acoustic attenuation 
r in metals. For a phonon with a wavenumber k and 
frequency u>k = sk, they are given by IJ 
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where 7 is the electronic scattering rate, and the di- 
mensionless coupling constant A is proportional to the 
electronic density of states Uq at the Fermi surface (for 
the isotropic case vq — m*pp /it 2 , m* is the effective 
electron mass) and to the squared deformation potential 
Cife- The deformation potential describes the change in 
the spectrum of electrons subject to lattice deformation 
e(p, r, t) = £ (p) + 0fc(p)Uifc(r, t), where u lk is the strain 
tensor. Equations (jl|) give the correct answers in various 
known regimes: for the sound attenuation in the hydro- 
dynamic limit (ujk "C 7 and fc — > 0), for the zero-sound 
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(u>k ^ 7 and k — > 0), for the Landau damping in the 
ballistic limit (kvp ^> \ojk + In the latter case, both 
the sound velocity shift and attenuation are small by the 
adiabatic parameter s/vp contrary to the results of the 
Frohlich model. Note also that Eqs. ([!]) show harden- 
ing of the phonon frequency due to the elect ron-phonon 
interaction in contradiction with Migdal's result. 

For the case of optical phonons, two known types of 
the electron-phonon interaction are the deformation po- 
tential and the interaction with the electrical polarization 
induced by optical vibrations. We consider here for sim- 
plicity a cubic crystal with two different atoms in a unit 
cell. Then, there are three optical modes, and the in- 
teraction with the induced polarization has the Frohlich 
form: 



e(p,r,t) =e (p) + C(p)V-u(r,t) 



(2) 



where the scalar function C(p) of the electron momentum 
is the coupling with the optical displacements u. In order 
to compare our results with previous ones, here we con- 
sider the interaction in the same form (||) as in Refs. 0, 
||. One can see that principal characteristic features of 
the phenomenon are retained for the deformation inter- 
action, 5e — £t(p)tti(r, t), where the coupling is a vector 
function. The distinction is that only the longitudinal 
mode interacts with electrons in the case of induced po- 
larization (Q), and the interaction approaches zero in the 
long-wave limit. Therefore, we concentrate on the propa- 
gation of the longitudinal mode along the symmetry axis 
when this mode is not mixed with transverse ones. Note 
also that the electric field plays an important role espe- 
cially when the different atoms are in the unit cell so 
that the dipole moment is excited under the atom vi- 
brations. At last, the optical phonons have always the 
so-called natural width r nat ~ vo^/m/M. The natural 
width results from decay processes into two (or more) 
acoustic or optical phonons, which are possible even at 
zero temperature. 

The main point of the theory is the equation of motion 
in the long-wave approximation (k <§C 1/a) for the Fourier 
components of the optical-phonon displacement Uj : 
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C(p)<S/ P (k,a,), 



where Ej is the electric field associated with vibrations, 
N is the number of unit cells in 1 cm 3 , M' is the reduced 
mass of two atoms in the unit cell, and Z is the effective 
ionic charge. The nonperturbed phonon frequency ojk 
should be considered in the absence of the electric field 
and without any nonadiabatic corrections. In the long- 
wave limit, we can roughly describe it as lo\ = uj 2 ± s 2 k 2 
with the magnitude of s being of the order of the typi- 
cal sound velocity in metals. The last term in Eq. (0) 



presents the driving force from the nonadiabatic elec- 
tron system due to the deviation 6f p (r, t) from the local- 
equilibrium distribution function f [e(p,r,t) — fx]. 
Then, we have the Boltzmann equation 

- i(w - k • v)Sf p (k, u ) + 7 [Sf p (k, uj) - (Sf p (k, W )>] (4) 



= -[wC(p)k-u(k,w) H-ev-E] 



de ' 



in the approximation of relaxation rate 7, which holds 
at low temperatures when the electron-impurity inter- 
action dominates as well as at temperatures higher then 
the Debye temperature when the phonon-phonon colli- 
sions can be considered as elastic. The term in the angle 
brackets in Eq. (Q), which denote the average over the 
Fermi surface, 
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arises from the out-term in the collision integral. No- 
tice, that the condition (C(p)) = should be fulfilled, 
because the number of electrons in the local-equilibrium 
state /o[e(p, r, t) — /x] is conserved. With the help of the 
Maxwell equations, the electric field is expressed in the 
terms of polarization P as follows: 
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-47rk(k-P)/A; 2 



(5) 



provided that phonons are excited in the optical region 
k ^> lo I 'c where the wavevector is determined by the in- 
cident light k ~ u;W/c and the frequency ui is of the 
order of the optical phonon frequency ujq. It is seen that 
the electric field is longitudinal and only the longitudi- 
nal component of polarization P z (k is taken along the 
z-axis) plays a role being related to the phonon displace- 
ment and the electric field by the equation 
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Sfpfow), (6) 



where the first term is caused by ionic motion, a is the 
polarizability of filled bands, and the last term is the car- 
rier contribution defined by the variation of the electron 
density, = ~~ik ■ . 

Equations (|)-(|) give the complete system of our 
problem. The Boltzmann equation (||) has the solution 
in the form 



-x P (k, w )-, 



where 



Xp (k, w) =i[ev-E + w((p)k • u + 7 (x P (k, w))]/A, (7) 
( Xp (k, w)) = i([ev ■ E + c.C(p)k • u]/A)/(l - i( 7 /A)), 

and we set A = u> — k • v + 27. 
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Using this solution, we obtain the polarization (g) and 
rewrite the electric field (||) in terms of the longitudinal 
displacement u z : 



£ e (k,u))E = -&nfifdu z , 



(8) 



where we introduce the field-displacement response func- 
tion 



O) 



The electron contribution into the dielectric function has 
the known form: 



4ire 2 is (v z /A) 
"fc(l-i( 7 /A})' 



(10) 



where the high-frequency permittivity Erx, = 1 + Aira. 

Now, we consider the equation of motion (||) using the 
solution of the Boltzmann equation (]?]). The term pro- 
portional to u z of the driving force can be included in 
the phonon frequency: 



Co 2 = uji - iujT nat + 



v uk 2 / 7 C 2 (P)\ , *7<C(P)/A> 5 



M'N y A 1 1 - 1(7/ A) 
so that Eq. @ reads 

(uj 2 - lu 2 )u z = ZE/M', 
where the renormalized ionic charge 

s 7 e^ofc ^ C(PK \ 1 ■ (^/A)(C(p)/A) ^ 



(11) 
(12) 



Using the condition (C(p)) = we obtain Z = fifd/N. 

Then, we can express the displacement u z from Eq. 
( |l2| ) in terms of E and, substituting into Eq. (||), obtain 
the dielectric function of the electron-ion system: 



e(k,Lu) =e*- e) (k,u) + 4nNZ' 2 /M'(u>' 2 -uo 2 ). 



(13) 



The frequency of the longitudinal mode, ui — wloj is 
defined by the condition e(k,tv) — 0, i.e., 



u 2 = 6j 2 + 4TrNZ 2 /M'e e (k, u). 



(14) 



In the absence of free electrons, the density of states 



= and Eq. rtTl) gives for the LO mode w£ 



LO 



tot 



iu)L T nat : where = AttNZ 2 /e^ is the squared 
ion-plasma frequency of the order of uJq. For the TO 
mode, when the electric field E = 0, we obtain w^o = 

lo 2 -^ fc r nat . 

Free electrons in metals make the large contribution 
into the dielectric function [see, Eq. (|l0|)]. Expanding in 
powers of k we have in the zero-order 



e^(k,Lu)-e c 



-SooUJ 2 e /uj(uj + 17), (15) 



which corresponds simply to the Drude conductivity with 
the electron-plasma frequency 



WpB 37T 2 £, 



vcISf- 



For large kvF > |w + the electron contribution ( |l0| ) 
describes the Debye screening: 

k 2 

e {e) {k,u) -600= e ooyJ(l + inuj/2kv F ), 

where a term of the order of w/kv is kept and the Debye 
parameter fcg = 47re 2 ^o/£co- 

Therefore, we can solve Eq. jl4| ) for w^o <C w pe , using 
the iteration procedure. To a first approximation, we 
have 



2 2 ■ T^r 



k 2 uj k v / C 2 (P) \ fpi^oo 
M'N ^ A ' e( e )(fc,w fe ) 



(16) 



In the case of small kvp < \u>h + vy\, expanding in 
fc-powers, we obtain the solution 

2 2 • F nat / 1 • \ UJ pi 1 AcJfcS 2 2 . , 

w LO= W fe _iW fc r -Wfc(w/c+«7)^-H ; — —k , (17) 

^pe Wfc + Z7 

where the dimensionless coupling constant A = 
{C, 2 {p))vq/ ps 2 contains the factor appm*/m and the 
metal density p. 

In the case of large k, expanding in \uok + i^/\/kvF, we 
obtain 



J LO 



.mVkS 2 k 
2v F 



A + 



(18) 

where the coupling constant A is defined, when 
the asymptotic value of the integral is calculated: 
^o(C 2 (p)/A)/Af'7V = -ins 2 X/2kv F . Note, that the 
value of A vanishes in the isotropic case due to the con- 
dition (C(p)} = 0. 

Hence, the squared frequency of the longitudinal opti- 
cal mode is essentially less (by the factor u 2 ^ than for 
insulators, due to screening of the electric field by free 
electrons. The additional phonon softening, width, and 
dispersion in Eq. ( |l7| ) involve the adiabatic parameter 
(wo/wpe) 2 ~ m/M, and they are small compared with 
r nat loq^ m/M . In the region, where Eq. (|l^) is valid, 
we see the additional fc-dependent width (terms in the 
parentheses), which is comparable with r nat . Here, the 
A-term conditioned by the elect ron-phonon interaction 
is similar to the damping of the acoustic phonons [see 
the first formula in Eqs. (|j])]. Now we omit the shift 
containing the small factor (s/vf) 2 - The second term 
in the parentheses as well as the last term in Eq. (p"8]), 
giving the fc— dispersion, are induced by screening and 
A— independent. Since w^/Zcq ~ s 2 , this dispersion has 
the typical value for the phonon branches. 
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Let us rewrite in our notations the respective results For the conducting phase, the longitudinal mode softens 
of Ref. S, Eqs. (9) and (7), retaining only main terms: and obtains the additional damping in comparison with 

the insulator state. Then, using Eq. (|l7|), we calculate 
2 _ 2 w j>« v^pi + w fc ) i / iq s ]l2t the collision rate 7 which is consistent with the value 

^o-^k- o (i_2Aln2)' ltv f<^, t iJ J obtained from the conductivity. 



(20) 



Comparing with Eqs. (|l7|)-(|l8|) we see that the fc- 
dispersion coincides practically. Next, we agree that 
the contribution vanishes from the frequency of the 
LO mode. The softening and damping due to both the 
electron-phonon scattering (7) and the phonon decay- 
processes (r nat ) were ignored in Ref. ||. Concerning 
the electron-phonon interaction A, the reason of dis- 
agreement was discussed in the beginning of the paper: 
this is shortcomings of diagram technique based on the 
Frohlich model. But the most essential difference is the 
/c-dependent width in Eq. (|l8|), which is missed in Eq. 

©■ 

In conclusion, let us make several remarks. The case 
of the large A:- values ([l8|) is most interesting because the 
electron-phonon and electrodynamic contributions into 
the phonon width (first and second terms in the parenthe- 
ses, respectively) can compete. The result depends on the 
Debye screening and the wavevector k. In Raman exper- 
iments, the parameter kvp/\ojQ + i"f\ ~ w^'Vf/ojJq ~ 0.3 
if ~ 10 4 K and lu ~ 10 2 K, and for metals v F ~ 10 8 
cm/s. Therefore, the high incident-light frequency or 
neutron experiments are desirable. It is more simply 
to observe the electrodynamic effect in semiconductors, 
where k$ is smaller. We have an example of such ex- 
periments in works jOJ], where the metal-insulator tran- 
sition was observed in the GaN crystal under pressure. 
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